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Abstract. We show that the number of cycles in a random per- 



(N 

o 

p^ ' mutation chosen according to generalized Ewens measure is nor- 

. , mally distributed and compute asymptotic estimates for the mean 

C$ ' and variance. 

Si 

1. Introduction 

Let S n denote the permutation group on n letters. For each permu- 
Q^ ■ tation a G S n , we write Cj(a) for the number of disjoint cycles of length 

j in a. For any permutation, we let K 0n (a) := Y^j=\ c j{ a ) denote the 
^ ', number of cycles in a. 

We are interested in the statistics of permutations produced in a ran- 
dom way. Random (uniform) permutations and their cycle structures 
have received much attention and have a long history (see e.g. the first 
chapter of [1] for a detailed account with references). The literature on 
the topic has grown quickly in recent years in relation to mathematical 
biology and theoretical physics, where models of non-uniform permuta- 
tions are considered (see e.g. [21 El HI IS])- We will restrict our attention 
to random permutations with cycle weights as considered in the recent 
work of Betz, Ueltschi and Velenik [I] and Ercolani and Uelstschi [5]. 
pq . These are families of probability measures on S n that are constant on 

conjugacy classes with the distribution 






Cj{a) 






3=1 



where 9j > is a given sequence of weights and h n is a normalization 
constant. If 9j = 1 for all j then this is the uniform measure on S n , 
while if 6j = 6q is constant then this gives Ewens measure, which plays 
an important role in mathematical biology. 

A situation of interest which appears in the study of the quantum 
gas in statistical mechanics is when the asymptotic behavior of 6j is 
fixed for large j (see [1] and [5]). Natural important historical questions 
arise on the behavior of Cj(a) or Ko n (a). For instance, it is known that 
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under the Ewens measure, or in special cases of weighted random per- 
mutations, the cycle counts (cj(a))j converge to independent Poisson 
distributions (see [TJ for the Ewens measure and [5] and [7J for weighted 
random permutations). The case of K 0n (a) is in fact more delicate and 
less results are available in the general case with cycle weights. It 
is well known that under the Ewens measure Ko n (a) satisfies a cen- 
tral limit theorem (see pQ for details and historical references). The 
methods used in this case are very probabilistic and rely on the Feller 
coupling. However, the Feller coupling does not exist in the model of 
random permutations with cycle weights. Ercolani and Ueltschi ([5]) 
used generating series and refined saddle point analysis to obtain some 
asymptotic estimates for the mean of K 0n (a) in some special cases but 
were not able to prove any central limit theorem. In [7J the second 
and third authors used generating series and singularity analysis to 
prove a central limit theorem and some large deviations estimates in 
the cases where the generating series exhibit some logarithmic singu- 
larities, but the important cases corresponding to subexponential and 
algebraic growth of the generating series were still open (see the corol- 
laries below for a more precise statement). In this paper we propose 
yet another, but more elementary, approach based on Cauchy's integral 
theorem for analytic functions to solve these problems. 

More precisely, with a sequence 9 = {9j} c *L 1 fixed, we write 

oo „ 

fc=i 
for the indicated generating function. We will always assume that 
the series for gg converges in a neighborhood of the origin. We will 
also require that gg satisfies a technical condition which we call log- 
admissibility, which will be defined in Section |2j 
Our main result is the following. 

Theorem 1.1 (Central Limit Theorem for K 0n ). Suppose that gg(t) is 
defined for t e [0, 1) with gg(t) — > oo as t — > 1. Suppose further that 
g is log- admissible (see Definition \2.1\) . Then there are sequences fi n 
and o~„ such that 

K{) n — fj, n 



' n 



0~n 

converges to a standard normal distribution. 

The mean and standard deviation can also be explicitly computed. 
For now we will state explicit results for two cases of interest, and defer 
the general result to Section |3j 

Corollary 1.2. Let g{t) = 7(1 — t)^ 13 for some (3 > and 7 > 0. Then 
K Qn converges asymptotically to a normal distribution with mean 

/Un = ( / 3 7 )-^+iV/(/3+i)(i + o( i)) 
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and variance 

°i = ((Pi)- 1 -w+ lnc/^w+w/^ci + o(i)). 

Corollary 1.3. Let g(t) = exp(l — t)~P for some (3 > 0. Then K 0n 
converges asymptotically to a normal distribution with mean 

n 

^ = (logn)i+i/g (1 + ° (1)) 

and variance 

ff " = (2 + rl) (lognV^ 

Our approach relies in the following well-known calculation of the 
moment generating function for Kq u . 

Proposition 1.4. We have the power series identity 

oo 

exp(e~ s g e (t)) = ^h n Eexp(-sK 0n )t n 

n=0 

for all sGK, with the conventions that ho = 1 and Koo = 0. 

This proposition follows immediately from Polya's enumeration the- 
orem with a small calculation. More details can be found for instance 
in pj Section 4] 

Remark. In [7j the characteristic function Kexp(itKo n ) of Ko n was con- 
sidered. In our new approach it is crucial to rather consider the Laplace 
transform Eexp(— sK 0n ) where the variable s is real in order to be able 
to evaluate the relevant contour integrals. 

The outline of this article is as follows. In Section [2] we define log- 
admissible g(t) and derive a formula for the coefficients of its generating 
function. In Section [3] we use the formula to compute asymptotics for 
Eexp(— sK 0n ), of which Theorem 11.11 is a direct consequence. We also 
prove Corollary 11.21 and Corollary 11.31 Finally, in Section H] we show 
how the proof of Theorem 11.11 can be modified to give large deviation 
estimates for K 0n . 

Notation. We will also freely employ asymptotic notation as follows; 
let /, g, and h be arbitrary functions of a parameter n. Then we write 
/ = 0(g) to indicate the existence of a constant C and threshold n 
such that for all n > n , \f(n)\ < C\g(n)\; the constant and threshold 
may be different in each use of the notation. We also write / = h+0(g) 
to indicate \f — h\ = 0(g). We similarly write / = o(g) to indicate 
that 

lim M = 0. 
n->oo g(n) 

It is convenient to also employ Vinogradov notation: we write / < g 
(and equivalently g > /) for / = 0(g). 
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2. Estimates for the moment generating function 

We shall now introduce the class of functions gg(t) where we can 
compute the asymptotic behavior of Eexp(— sK Qn ). 

Definition 2.1. Let g(t) = Yl^=o9nt n be given with radius of conver- 
gence p > and g n > 0. We say that g(t) is log-admissible if there 
exist functions a,b } 5 : [0, p) — > R + and R : [0, p) x (— n/2, n/2) —> W + 
with the following properties. 

approximation: For all \ip\ < 5(r) we have the expansion 

2 

g(re ilp ) = g(r) + i<pa(r) - — b(r) + R(r, <p) (2.1) 

where R(r,(p) = o((p 3 5(r)~ 3 ) and the implied constant is uni- 
form. 

divergence: a(r) — > oo ; b(r) — > oo and 5(r) — > as r — > p. 

width of convergence: For all e > 0, we have e5 2 (r)b(r) — 
logfe(r) — > oo as r — > p. 

monotonicity: For all \tp\ > S(r), we have 

n g {re iv ) < ^g{re ±iS{r) ). (2.2) 

These properties can be interpreted as a logarithmic analogue of 
Hayman-admissibility [B] Chapter VIII. 5]. In fact, if g(t) is log- admissible 
then expg(t) is Hayman-admissible. We have also introduced the e 
term in the width condition, which is required for uniformity in the 
error term of Proposition 12.21 

The approximation condition allows us to compute the functions a 
and b exactly. We have 

a(r)=rg'(r), (2.3) 

b(r) =rg'(r)+r 2 g"(r). (2.4) 

Clearly a and b are strictly increasing real analytic functions in [0, p). 
The error in the approximation can similarly be bounded, so that 

R( r > <P) ^ r g'( r ) + 3r 2 g"(r) + r 3 g'"(r). 

Note that for Ewens measure (6j = 1 for all j > 0), we have g(r) = 
— log(l — r), so we can compute 



^( r ) = Ti v? an d -^( r 5 V) 



qr> ■" I iyi 



(l_ r )2 ^^> (l- r ) 3 

Therefore g(r) = — log(l — r) is not log-admissible and we cannot apply 
this method to such distributions. 

We will frequently require the inverse function of a on the interval 
[0,p), and define r x to be the (unique) solution to a{r) = x there. It is 
easy to see that r x is strictly increasing real analytic function in x and 
r x — > p if x tends to infinity. 
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We now define for s G K the sequence G n)S by 

oo 

J2G n , s t n = exp(e- s g(t)). (2.5) 

n=0 

If g(t) is log-admissible, then we can compute the asymptotic behav- 
ior of G ntS for n — > oo. The generating function in Proposition 11.41 has 
the same form as (12.51) and we thus can compute also the asymptotic be- 
havior of the moment generating function for K§ n if gg is log-admissible. 

Proposition 2.2. Let s6R and let g be log-admissible with associated 
functions a, b. Let further be r x to be the (unique) solution of a{r) = x. 
Then G nyS has the asymptotic expansion 

G n , s = 4=e s/2 r e -46(r eSn )" 1 / 2 exp( e -^(r e , n ))(l + o(l)) 



as n —¥ oo with the implied constant uniform in s for s bounded. 

Proof. We apply Cauchy's integral formula to the circle 7 centered at 
of radius r. We get 

1 f ,„,,., dz 



G n>s = — / exp(e s g{z)) 

27U J 



7 Z 



exp(e s g(re %tp ) — imp) dip 



2vrr n J _ r - 

with r = r e s n . We now split the integral into the parts [— 8(r),8(r)] 
and [—71, — S(r))U(S(r),ir]. We first look at the minor arc [— <5(r), S(r)]. 
By hypothesis on g we can expand in (p, giving 

/S(r) 
exp(e~ s g(re l<p ) — imp) dtp 
■S(r) 

exp(e~ s (g(r) + iipa(r) - b(r)ip 2 /2 + o(ip 3 5(r)~ 3 ) - imp)) dip. 

-S(r) 

We have e~ s a(r e s n ) = n since r = r e s n , which cancels the linear terms 
in (p. We get 

h = [ exp(e" s (^(r) - b(r) V 2 /2 + <p 3 o(5(r)- 3 ))) dip. 

J-8(r) 

We now observe that (p 3 o(5(r)~ 3 ) = o(l) on [S(r), S(r)] as r = r e s n —$■ 
p, or equivalently as n — > 00. Rearranging, we get 

h = exp{e~ s g(r)) / exp(- e - s 6(r)^ 2 / 2 ) d<p{l + o(l)). 

J-5{r) 

By assumption on the width of convergence of g, the integral converges 
to 

/•<5(r)e- s / 2 b(r) 1 /2 

e s ' 2 b(r)- 1/2 / exp(-x 2 /2) dx = v / 2Te s/2 6(r)- 1 / 2 (l + o(l)) 

J -S(r)e- s / 2 b(r)V 2 
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SO that 

h = V2^exp(e- s g(r))e s/2 b(r)- 1/2 (l + o(l)). 

Next we estimate the integral over the major arc. By the monotonic- 
ity assumption on g, we have 

J 2 := I exp{e' s g{re i(p )-inip) dip < 27rexp(3?(e~ s #(re^ (r) ))) 

J[-ir,-5(r)]U[5(r),n] 

We can apply the approximation for g at ip = S(r), giving 

n g (re i5 ^) = g(r) - 5 -^b(r) + o(l) 
Collecting terms and rearranging, 

I 2 < 2tt exp(e~ s g(r))b(r)- 1/2 exp(-e- s 5(r) 2 b(r) /2 + -logfe(r)) 

= o(exp(e- s g(r))b(r)- 1/2 ) (2.6) 

where at the last step we used the width of approximation for g. 
Combining the estimates for I\ and I 2 , we find that 

G n , s = -^=r- n n e s/2 exp{e- s g(r e s n )){l + o{l)). 
V In 

where the error term is uniform in s for s in a fixed compact set. 
Note that all error-terms in this proof are uniform in s for s in a fixed 
compact set. □ 

Note that the e in the definition of log-admissibility is required to 
make the error in (J2.6P uniform in s. 

3. The total number of cycles 

We are now ready to compute the asymptotic number of cycles in a 
random permutation as described in the introduction. We will restrict 
our attention to those examples in [5] where the limiting behavior was 
not known, namely where the generating function g$ is of the form 

g e (r) =7(1 - r)^ 

or 

g e (r) = exp(l - r) _/3 . 

We will refer to such functions as exhibiting algebraic and sub-exponential 
growth, respectively. 

We begin by observing that a formula for the moment generating 
function of K 0n can be determined from Proposition 12.21 

Corollary 3.1. Let s G M, ge(t) be log-admissible with associated func- 
tions a, b. Let further r x be the (unique) solution to a(r) = x. 

Then 

1 

hn= 7tr=rr=^ ex P(^(r n ))(l + o(l)) 

A/27r6(r n )r^ 
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and 

Eex (sK ) = c s/2 ( Tn X ex P( e ~^( r ^ s ")) ( fe ( r ") V /2 q | (-y\\ 

\r e s n ) exp(g e (r n )) \b(r e s n ) ) 

where the error terms are uniform in s for s in a fixed compact set. 

Proof. By Proposition II. 44 we have 

h n Eexp(-sK 0n ) = G n>s . 

Apply Proposition 12.21 at to get the desired formula for h n , and apply 
it again at s to find the formula for Eexp(— sKo n ). D 

3.1. A simple example. Before we consider more complicated func- 
tions, we will illustrate the method with gg given by the equation 

1 



99(f) 



l-t 



This generating function corresponds to the sequence 9k = k. 

Our first step is to compute the moment generating function of Kq h 
by finding an asymptotic expansion of the formula in Corollary 13.11 

Proposition 3.2. Let gg(t) = (1 — t)~ l . Then gg is admissible, 
h n = ^n- 3/4 exp(2 v / ^(l + o(l))), 

V47T 



and 

Eexp(-sK 0n ) = e- s/4 exp(2 v / ^(e" s/2 - l)(l + o(l))) 
where the errors are uniform in s for s in a fixed compact set. 

We will prove Proposition 13.21 in a moment. We first see how to use 
this result to prove a central limit theorem for K 0n . 

Corollary 3.3. Let g e (t) = (1 - i)" 1 . Then 

K 0n -n l l 2 d 

2-l/2 n l/4 

where N is the standard normal distribution. 

Proof of Corollaru \3.3[ It suffices to show that the moment generating 
function of the renormalized cycle count converges to e s I 2 for bounded 
s (in fact, we only need this for s sufficiently small, but our theorem 
gives a stronger result). We apply Theorem 13.21 at s/(2~ 1//2 n 1 / 4 ) to find 

E exp(-* 2 -i^i/ 4 ) = ex P (-V2sn 1 ^ + s 2 /2 + O^n^)). 

Here we used the uniformity of the error for bounded s. Multiplying 
both sides by exp(-y/2sn 1 / 4 ) and letting n tend to oo completes the 
proof. □ 
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Proof of Theorem \3.B. We first show that gg(t) = (1 — t)" 1 is admissi- 
ble. The monotonicity condition is obvious. We compute 

a(t) = a^F 

and 

t It 2 

m = (I" If + (137)3 

and note that a, b — > oo as t — > 1 from the left. It suffices to choose 
a function 5(t) that satisfies the remaining hypotheses. We observe 
that the width condition on 5 is satisfied if S(t) < (1 — i) 3 / 2 ^ for some 
7] > 0. Likewise, for the error in the approximation condition to be 
satisfied we need 5(t) 3 (l — t)~ A = o(l). It therefore suffices to choose 
S(t) = (1 - t) a for any 4/3 < a < 3/2. 

We first calculate r x . By definition we have 

r x 
(l-r,) 2 " X ' 

which can be inverted to find 

r x = l-x- l ' 2 {l + o{l)). 



We then compute 
and 



g(r x ) = yfc(l + o(l)) 



b(r x ) = 2x 3 / 2 (l + o(l)). 

With the approximation (1 — r]) n ~ exp(— r]n) for rj sufficiently small, 
we apply Corollary 13.11 to find 

h n = -^n- 3 / 4 exp(2v^(l + o(l))) 

V47T 



and 

Eexp(-s J fs:on) = e- s/4 exp(2v^(e- s/2 - 1)(1 + o(l))) 
as required. □ 

3.2. The general case. The previous calculation suggests how to 
transform the formula in Corollary 13.11 into a form that is easier to 
manage. We will restrict our attention to those functions g where the 
induced functions r x satisfy a family of inequalities 

|r«x fe |<|rrV^)| (3.1) 

for all x sufficiently large and 2 < k < 4. This is easy to verify in 
practice and avoids some technical details for functions g which diverge 
slowly at 1 (i.e. slower than (1 — r)~ e for any e > 0). 
It is convenient to define the functions 

rik(x) := (-l) fc 



d k s 



log(rw) 
=o 
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For example, this gives 

/ / // / / \ 2 

m(x) = —x— and r] 2 (x) = -r]i(x) + x 2 I — + ' ' 
r x yr^. 

Proposition 3.4. Let g be log- admissible and r x be defined as above, 
satisfying condition \3. 11 Fix M > 0. Then for every —M < s < M, 
we have 

\ogEexp(-sK 0n ) =-g(r n )(l + o(l))s 

+ {g{r n )+nr ]l {n)){l + o{l))s 2 /2 

where 

£(n) = sup g(r x ) +x\r] 1 (x)\. 

x£[e~ M n,e M n] 



Proof. We use Corollary 13.11 and expand each factor with respect to s. 
We start with the first term. Taking logarithms, we find 

log(r n r~} n ) n = nlog(r n ) - nlog(7>„) 

= nr]i(n)s — nrj 2 (n)s 2 /2 + 0(^i(n)s 3 ). 

where 

£i(n) := sup |%(ar)|. 

x£[e~ M n,e M n] 

Next we consider exp(e~ s g(r e s n ) — g(r n )). We use x = a(r x ) = r x g'{r x ) 
and get 

o r 

j^-g(r e s n ) = g'{r e s n )r' ean e s n = e 2s n 2 ^^ = -e s n ■ r) X (e s n). 

OS T e a n 

This gives 

e~ s g(r e s n ) - g(r n ) = - (71771(71) + g{r n ))s 

+ (71771(71) + ni] 2 (n) + g{r n ))s 2 /2 
+ 0(Un)s 3 ) 
where 

£2(71) = sup g r (x) + x\i]i(x)\ + x\rj 2 (x)\ + x 1 773(2;) I . 

x£[e — M n,e M n] 

Finally we consider {b{r n )b{r e s n )~ 1 ) 1 ^ 2 . Writing b in terms of a, we 
get 

b{r x ) = a(r x ) + r 2 x g"(r x ) = x + r 2 x g"(r x ) 
Differentiating the defining equation x = r x g'(r x ) in x, we find that 

1 = r' x g'(r x ) + r x r' x g"(r x ) 
so that 

r x9 (>"*) = ^j-x 



10 

and thus 
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X 



-X 



Vi{x) 



log 



This then gives 
b{r n ) \ 1/2 



b{r e 



s/2- 



2r)x(n) 



r} 3 {n) 



r]l(n) 



2r)i(n) 2r] 2 (n) 



■+0(Un)s 3 ) 



where 



Un) 



sup 

x£[n,e s n] 



77 4 (Xj 



Tll(X) 



+ 



rj 2 {x)r} 3 (x) 



7/2(2) 



+ 



•ql(x) 



rif(x) 



We now apply the technical assumption to see that f]k{n) 



< 



Vi\n) 



for k = 2,3,4. In particular, we see that the (6(r n )6(r eSn ) -1 ) 1 / 2 is 
dominated by the other terms. Collapsing other redundant terms, we 
find that 



logEexp(— sK { 



on] = - g{r n )(l + o(l))s 

+ {g{r n )+n Vl {n)){l + o(l))s 2 /2 
+ 0(s 3 sup g(r x ) + x\i]i(x)\) 

xG [e~ M n,e M n] 

as desired. 

As above, this has the following immediate corollary. 



U 



Corollary 3.5. Let 9 be the defining sequence for a generalized Ewens 
measure. Suppose that gg is log- admissible andr x satisfies the technical 
condition IJ.il Then there are functions \x n and a n such that 

Kon 



/-'-,, 



On 



N(0, 1) 



where fi n , a n satisfy the asymptotics 

Hn = g{r n ){l + o(l)) 
and 



o- 



(g(r n ) + n Vl (n))(l + o(l)). 



Proof. The only thing to check is whether the coefficient of s 3 is bounded 
by (g( r x) + x\gi(x)\) 3 ^ 2 , but this is obvious. □ 

3.3. Computing g(r n ) and //i(n). We now give two examples of how 
to apply Corollary 13.51 to compute explicit asymptotics for given gg. 
First we prove Corollary 11.21 

Proposition 3.6. Let g$(t) = 7(1 — t)^ 13 for some (3 > and 7 > 0. 
Then gg is log-admissible, r x satisfies the technical condition \3 . 1\ and 
there are asymptotic expansions 



*w = MTr(i + o(i)) 
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and 

nT]l (n) = - [P7 > - n^ +l \l + o(l)). 

Proof. Admissibility follows once we construct an explicit 5(t). It suf- 
fices to find a function that satisfies the inequalities 

e5(t) 2 b(t)-\ogb(t) ^oo 

for all e > 0, and 

6(t) 3 R(t,<p) ->0. 
For t — > 1, we have the lower bound &(£) > 7ij g (1 — £) _/3 ~ 2 and the upper 
bound R(t,ip) < 7ii g (1 — t) _/3_3 . Thus we see that any 5 of the form 
S(t) = (1 - t)° with l + |<a<l + | suffices. 

We compute r n by inverting n = a(r n ) = f3^r n {l — r^ 13 " 1 , so that 

r n = l-(/3 7 n- 1 )^i(l + o(l)). 

The derivatives of r n can be approximated in an analogous way, so that 

lr (fc) l = Cr n^~ k 

I'n I v ^p,7"' 

and the technical condition is clear. We then estimate 

^(r„) = (n/3- 1 7 - 1 )^ T (l + o(l)). 
and 

(o \ 1/03+1) 

nm (n) = ~ 1 ! ^ /(/9+1) (l + o(l)) 
with our estimate for r n and r^. D 



' (i + ( io gn) -v S + (i + r')^^(i + 0(D)) 



Next we prove Corollary 11.31 

Proposition 3.7. Let ^(r) = exp(l — r) _/3 /or some (3 > 0. T/ien o^ 
zs log-admissible, r x satisfies the technical condition \3.1\ and there are 
asymptotic expansions 

9e{Tn> (lognY+W^ ' v ~ 6 " 7 ' v- ' '' , logx 

and 

nm(n) = - " (1 + (\ogn)-^ - (1 + /^^(l + o(l))) 

[iog n) i+i ' p logx 

Proof. First we verify that go is admissible. Monotonicity is obvious. 
We compute 

a(r) = rg'ir) = r/3(l — r) - ^" 1 exp(l — r)~^ 

and 

b(r) = rg\r) + r 2 g"(r) 

= r 2 (/3(/3 + 1)(1 - r)-^ 2 + /3 2 (1 - r)" 2 ^" 2 ) exp(l - r)^ 

= r 2 /3 2 (l - r)- 2/3 - 2 exp(l - r) _/9 (l + 0/j(l - rf). 
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These diverge at 1. Once we estimate 

R(r, ip) = rg'(r) + 3r g"(r) + r g"'(r) 

= r 3 /3 3 (l - r)- 3/3 " 3 exp(l - r)^(l + O p {l - rf) 

we see that it remains to choose any S that satisfies the pair of inequal- 
ities 

e5(r) V/3 2 (l - r)- 2fi - 2 exp(l - r)~ fi -> oo 

and 

r 3 /3 3 (l - r)- 3/3_3 exp(l - r) _/3 < <5(r)- 3 . 

Any 5 of the form 

5(r) = exp(a(l-t)- fj ) 

with 1/3 < a < 1/2 suffices. 

We next need an asymptotic approximation for r x . This is provided 
by the following lemma. 

Lemma 3.8. Let f(x) := (1 — r^.)^ 13 . Then we have the asymptotic 
expansion 

f(x) = logx — (1 + /3 _1 ) log log x — log/3 

+ ((logx)- 1 ^ + (i + /r 1 )^^! + o(i)) 

logx 

Furthermore, we have the estimates 



a; fc log a; ' x fc (log:r) 2 

Proof. Once we make the substitution /(x) = (1— r x ) _/3 in the equation 
a(r x ) = x, we see that / is implicitly defined by the equation 

x = 0(l-f(x)- 1 ^)f(x) 1+1 ^expf(x). 

We then substitute f(x) = logo; — (1 + /3 _1 ) log log x + w and observe 
that w = ((log*)- 1 //? + (1 + ^)*&HL)(l + (1)). 

For the estimates on the derivatives of /, we differentiate the defining 
equation for / to find 

1 = {f{xf^ + (3(f(x)^ - 1) + 0(f(x)W - l)f(x))f'(x) exp f(x). 

We can use the defining equation again to eliminate the exponential 
term, which gives us 

fix) = Ui - p-'d-m-^r' + i 

Jy) x y p-i(i-f( x )-i/p)-i + i + f( x y- 

This gives us the lemma for k — 1. For the higher derivatives, we 
differentiate by parts and apply our earlier asymptotics. □ 
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Note that this lemma also shows that r x satisfies the technical con- 
dition. We apply this formula to g(r n ) to get 

n n , n \-i/B , /i , a -i^°S^ogx i 



9\r r , 



1 + (logn)- 1 ^ + (1 + /r 1 )^^! + (1))) 



(logn) 1+1 /^ v logx 

and to nrji(n) to get 

nvi(n) = - " (l + ilognyW-il + p- 1 ) r^- (l + o(l))). □ 

{[ogn) l+1 ' p log a; 

Other g{t) can be computed in similar ways. Note that in the proof 
of Corollary II. 3| it was crucial to develop g(r n ) and nrji(n) beyond the 
first term; this reflects the reduced variance of the number of cycles 
when there most of the cycles are of logarithmic length. 

4. Large deviation estimates 

The method developed in the previous two sections actually gives 
more information than a central limit theorem. In fact, it was enough 
for us to show that 

E^^ = expAl + o(l))) 
a n 2 

for s arbitrarily close to 0, but our method applied to all s in a fixed 
compact set. In this section we will briefly indicate how to use this 
extra information to prove large deviation estimates for Kq u . 

Let M(s) denote the moment generation function for the renormal- 
ized cycle count; i.e. 

M(s) = Eex P (s K ° n ~^ n ) 

On 

and let A(s) = log M{s) denote its logarithm. We restate the corollary 
of Proposition 13.41 as follows. 

Proposition 4.1. There are functions o~\, £(n) such that the for all 
s = 0(a n ), we have the estimate 

A(s) = s 2 /2 + 0(an)cr~ 3 )s 3 . 

As an immediate consequence, we also have 

A'(s) = s + 0(an)a- 3 )s" 

and 

A"( S ) = l + 0(((nKl, 

Furthermore, A'(s) is monotone increasing (hence injective) for such s. 
Theorem 4.2. For all a = O(o~ n ) we have 

Kon-fMn _ < _ ^ + exp( _ fl 2 /2 + q^ + £fl)) 
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where 

5 = 0{i{n)a-*a) 

and the errors are absolute. 

Proof. Let X n := (K 0n — fi n )/cr n and let i] denote the pdf for X n . We 
define a pdf v s depending on s G 1 by 

1 

Then if Y is a random variable with pdf i/ s , we see that 

P(|X n - a\ < e) = M(s)Ee- sY l lY _ al<e 
We want to choose s so that Y has mean a. In fact, 

EY = M(s)- 1 EX„exp(sX„)] = A'(s) 

Therefore, because A' is injective we solve s = a + 0(£(n)<7~ 3 a 2 ). On 
the event that \Y — a\ < e, we see that e~ sY = e - sa +°( se ) so that 

F(\X n - a\ < e) = exp(-a 2 /2 + 0{aa 3 + ea))F(\Y - a\ < e). 

It is not hard to show that for s chosen so that a = EY, 

E\Y-a\ 2 = A"(s) 

so that by the second moment method, 

F(\Y - a\ < e) = 1 - e - 2 (l + 0(£(ny r ; 3 a)), 

and the result follows. □ 



References 

[1] Arratia, R., Barbour, A., and Tavare, S. Logarithmic combinatorial 
structures: a probabilistic approach. EMS Monographs in Mathematics. Eu- 
ropean Mathematical Society (EMS), Zurich, 2003. 

[2] Betz, V., AND Ueltschi, D. Spatial random permutations and infinite cycles. 
Comm. Math. Phys. 285, 2 (2009), 469-501. 

[3] Betz, V., and Ueltschi, D. Spatial random permutations and poisson- 
dirichlet law of cycle lengths. Preprint, 2011. 

[4] Betz, V., Ueltschi, D., and Velenik, Y. Random permutations with cycle 
weights. Ann. Appl. Probab. 21, 1 (2011), 312-331. 

[5] Ercolani, N., and Ueltschi, D. Cycle structure of random permutations 
with cycle weights, preprint, 2011. 

[6] Flajolet, P., and Sedgewick, R. Analytic Combinatorics. Cambridge Uni- 
versity Press, New York, NY, USA, 2009. 

[7] Nikeghbali, A., and Zeindler, D. The generalized weighted probability 
measure on the symmetric group and the asymptotic behaviour of the cycles. 
To appear in Annales de L'Institut Poincar, 2011. 



ON THE NUMBER OF CYCLES IN A RANDOM PERMUTATION 15 

Institut fur Mathematik, Universitat Zurich, Winterthurerstrasse 
190, 8057-Zurich, Switzerland 

E-mail address: kenneth.maples@math.uzh.ch 

Institut fur Mathematik, Universitat Zurich, Winterthurerstrasse 
190, 8057-Zurich, Switzerland 

E-mail address: ashkan . nikeghbaliOmath . uzh . ch 

SONDERFORSCHUNGSBEREICH 701, FAKULTAT FUR MATHEMATIK, UNIVER- 
SITAT Bielefeld, Postfach 10 01 31, 33501 Bielefeld, Germany 
E-mail address: zeindler@math.uni-bielefeld.de 



